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The transition process induced by the interaction of acoustic disturbances in the free- 
stream with boundary layers over a 5-degree straight cone and a wedge with blunt tips is 
numerically investigated at a free-stream Mach number of 6.0. To compute the shock and 
the interaction of shock with the instability waves the Navier-Stokes equations are solved in 
axisymmetric coordinates. The governing equations are solved using the 5 th -order accurate 
weighted essentially non-oscillatory (WENO) scheme for space discretization and using 
third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. 
After the mean flow field is computed, acoustic disturbances are introduced at the outer 
boundary of the computational domain and unsteady simulations are performed. Generation 
and evolution of instability waves and the receptivity of boundary layer to slow and fast 
acoustic waves are investigated. The mean flow data are compared with the experimental 
results. The results show that the instability waves are generated near the leading edge and 
the non-parallel effects are stronger near the nose region for the flow over the cone than that 
over a wedge. It is also found that the boundary layer is much more receptive to slow 
acoustic wave (by almost a factor of 67) as compared to the fast wave. 


I. Introduction 

I ransition from laminar to turbulent state in shear flows occurs due to evolution and interaction of different 
^ disturbances inside the shear layer. Though there are several mechanisms and routes to go from a laminar to a 
turbulent state, most of them are generally follow these fundamental processes; a) receptivity b) linear instability c) 
nonlinear instability and saturation d) secondary instability and breakdown to turbulence. 

In the receptivity process, unsteady disturbances in the environment such as acoustic waves and turbulence 
interact with the inhomogeneities in the geometry such as roughness and generate instability waves inside the shear 
layer. In quiet environments, the initial amplitudes of these instability waves are small compared to any 
characteristic velocity and length scales in the flow. In the linear instability stage, the amplitudes of these instability 
waves grow exponentially downstream and this process is governed by the linearized Navier-Stokes equations. 
Further downstream, the amplitudes of the disturbances become large and the nonlinear effects inhibit the 
exponential growth and the amplitude of the waves eventually saturate. In the next stage, these finite amplitude 
saturated disturbances become unstable and two- and/or three dimensional disturbances. This is called the secondary 
instability and beyond this stage the spectrum broadens, due to complex interactions and further instabilities, and the 
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flow becomes turbulent in a short distance downstream. In this paper we are concerned with the receptivity to slow 
and fast acoustic waves and the linear instability of hypersonic boundary layers over a cone and a wedge. 

II. Literature Review 

There have been a number of investigations conducted on the interaction of acoustic waves with supersonic 
boundary layers. The interactions of acoustic waves with a supersonic boundary layer using inhomogeneous stability 
equations were investigated by Mack 1 and Gaponov 2 . One important finding was that due to the interaction, the 
acoustic waves excite disturbances inside boundary-layer, which are much larger than that in the free stream. The 
interaction of stream acoustic waves with a non-parallel boundary layer was studied by Gaponov and Smorodsky 3 . 
The analysis and the calculations showed that the disturbances inside the boundary layer reach values significantly 
higher compared to that in the free-stream. It was also observed that there exists a critical Reynolds number where 
this excitation is the highest. 

For a supersonic boundary layer with sufficiently high Mach number to allow both first and second Mack 
modes 4 , Fedorov and Khoklov 5 considered boundary layer response to both the fast and slow acoustic waves. The 
boundary layer modes excited near the leading edge by the two acoustic waves can be referred to as Mode F and 
Mode S, for convenience. The work of Fedorov and Khoklov 6 and Fedorov 7 identified two receptivity mechanisms 
in this Mach number regime: (1) leading-edge receptivity and (2) inter-modal exchange between Mode F and Mode 
S. For adiabatic wall, Federov 7 found that receptivity to slow acoustic waves could be as much as 50 times the 
receptivity via the fast acoustic waves. Thus, the leading edge receptivity via the slow mode excitation is much 
stronger than in the case of inter-modal exchange. According to Fedorov 7 , this receptivity mechanism may gain 
significance in the highly cooled boundary layers. 

Ma and Zhong 8, 9 ’ 10 performed direct numerical simulation for a Mach 4.5 flat-plate boundary layer to 
investigate receptivity to fast and slow acoustic waves and the mechanisms of inter-modal exchange. They 
employed fifth-order accurate shock fitting method to solve the governing equations. Egorov, Fedorov and 
Soudakov 11 investigated a similar problem at a Mach number of 6.0 whereby simulation of receptivity to slow and 
fast acoustic waves and the effect of incidence angle on the receptivity were studied. 

The transition process induced by the interaction of acoustic disturbances in the free stream is numerically 
investigated for a boundary layer over a flat plate with a blunted leading edge at a free stream Mach number of 3.5 
by Balakumar 12 . The governing equations are solved using 5 th -order accurate weighted essentially non-oscillatory 
(WENO) scheme for space discretization and 3 rd -order TVD Runge-Kutta scheme for time integration. 
Balakumar 13 also investigated the receptivity of boundary layers over blunt flat plates and wedges at a free stream 
Mach number of 3.5 and at a high Reynolds number of 10 6 /inch. The linear stability result of his work showed that 
the bluntness has a strong stabilizing effect on the stability of two dimensional boundary layers. It was also revealed 
that the boundary layers on blunt wedges are far more stable than on blunt flat plates. Malik and Balakumar 14 
investigated the receptivity of supersonic boundary layers to acoustic disturbances at a free stream Mach number 
4.5. The results showed that the instability waves are generated near the leading edge region and that the boundary 
layer is much more receptive to slow acoustic waves by almost a factor of 20 compared to fast acoustic waves. The 
effect of the acoustic wave incidence angle was also investigated and it was found that the receptivity of the 
boundary layer on the wind ward side decreases when the incidence angle is increased. 

An experimental investigation was conducted on a 5 -degree half-angle cone in a conventional Mach 6 wind 
tunnel by Horvath et al 15 to examine the effects of facility noise on boundary layer transition. They checked the 
influence of tunnel noise on the transition onset points by comparing transition locations determined from their test 
to that previously obtained in a Mach 6 low disturbance quiet tunnel. 

Here, we employ a fifth order weighted essentially non-oscillatory (WENO) scheme for spatial discretization and 
use third order total variation diminishing (TVD) Runge-Kutta scheme for time integration to solve for the 
hypersonic boundary layer receptivity problem. For this study, we have selected the flow conditions and geometry 
used by Horvath et al 15 . Cone has a small nose radius (R n =0.001-in) and the flow around the leading edge is resolved 
by using a sufficiently dense grid. We assume adiabatic wall conditions for steady flow computations. The 
objectives are to understand the receptivity process near the leading edge of a cone and to estimate the receptivity 
coefficient of the instability waves generated near the leading edge. Computations are performed to determine 
whether the slow or the fast acoustic waves are more efficient in generating the instability waves. Also comparisons 
of computed shock stand off distances with the experimental results are shown. To compare the receptivity process 
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between the axisymmetric and two-dimensional geometries computations are performed for the hypersonic flows 
over a cone and a wedge. 


III. Governing Equations 

The equations solved are the conservative unsteady compressible two-dimensional axisymmetric Navier- Stokes 
equations 

dQ dF dG dF v dG v 0 
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where Q is the solution flow field vector, F and G are the axial- and radial-direction inviscid flux vectors given by 
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and F v and G v are the axial- and radial-direction viscous and heat conduction flux vectors given by 
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The vector S contains viscous flux and source terms associated with the axisymmetric geometry. In 2D 
formulation, S will be zero and shear stresses and heat fluxes are 
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In axisymmetric formulation source term, shear stresses and heat fluxes have the following form 
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There exist a singularity along the axis of symmetry y=0. After applying L’Hopital rule and taking limit as y 

dll dp dT dp 

goes to zero using the symmetry conditions = = = = 0, V— 0 at y=0 , source term along the axis 
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Applying same procedure to shear stresses and heat fluxes we get new expressions for the shear stresses at the 
axis y=0. 
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Here (x, y) are the two-dimensional and (x, y, 0) are the axisymmetric coordinates and (u, v) are the 
corresponding velocity components, p is the density,^ is the pressure. E is the total energy given by 


E = e + 
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Here e is the molecular internal energy and T is the temperature. 

The viscosity (ju) is computed using Sutherland’s law and the coefficient of conductivity (k) is given in terms of 
Prandtl number (Pr). The variables p, p, T and velocity are non-dimensionalised by their corresponding reference 

variables p^, p x , and V RT respectively. The reference value for lengths is computed by -Jv^Xq / U ^ , where x 0 

is the location of the beginning of the computational domain in the stream wise direction. For the computation, the 
equations are transformed from physical coordinate system (x, y) to the computational curvilinear coordinate system 
(§,r|) in a conservative manner and the governing equations become 

dQ dF dG dF v dG v - 

dt d% df] d § df] 


The components of the flux in the computational domain related to the flux in the physical domain by 
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IV. Solution Algorithm 

The governing equations are solved using 5 th order accurate weighted essentially non-oscillatory (WENO) 
scheme for space discretization and using third order total-variation-diminishing (TVD) Runge-Kutta scheme for 
time integration. These methods are suitable in flows with discontinuities or high gradient regions. These schemes 
solve the governing equations discretely in a uniform structured computational domain in which flow properties are 
known at the grid nodes. WENO scheme approximate the spatial derivatives in a given direction to a higher order at 
the nodes, using the neighboring nodal values in that direction. TVD-RK scheme integrates the resulting equations 
in time to get the point values as a function of time. Since the spatial derivatives are independent of the coordinate 
directions, the method can easily add other dimensions. It is well known that approximating a discontinuous 
function by a higher order (two or more) polynomial generally introduces oscillatory behavior near the 
discontinuity, and this oscillation increases with the order of the approximation. The essentially non oscillatory 
(ENO) and the improvement of these WENO methods are developed to keep the higher order approximations in the 
smooth regions and to eliminate or suppress the oscillatory behavior near the discontinuities. They are achieved by 
systematically adopting or selecting the stencils based on the smoothness of the function, which is being 
approximated. Shu 16 explains the WENO and the TVD methods and the formulas. Atkins 17 gives the application of 
ENO method to the N-S equations. Balakumar 18 describes in detail the solution method implemented in this 
computation. 

At the outflow boundary, extrapolation boundary condition is used. At the wall, the simulation uses viscous 
conditions for the velocities and an adiabatic condition for the temperature. The density is computed from the 
continuity equation. In the mean flow computations, the simulation prescribes the free-stream values at the outer 
boundary, which lies outside the bow shock. The blunt cone is assumed to align with the free stream. In the unsteady 
computations, it superimposes the acoustic perturbations to the uniform mean flow at the outer boundary. The 
procedure is to first compute the steady mean flow by performing unsteady computations using a variable time step 
until the maximum residual reaches a small value -10" 11 . These computations use a CFL number of 0.2. The next 
step is to introduce unsteady disturbances at the outer boundary of the computational domain and to perform time 
accurate computations to investigate the interaction and evolution of these disturbances to downstream. Since we 
use very fine spatial grid to resolve the leading edge region, these computations require very small time step which 
is taken as the minimum time step allowable for a CFL number of 0.2. 

The grid is generated using analytical formulae. The grid stretches in the r| direction close to the wall and is 
uniform outside of the boundary layer. In the § direction, the grid is symmetric about the leading edge and very fine 
near the nose and is uniform in the flat region. The outer boundary that lies outside the shock follows a parabola so 
the boundary layer growth could be captured accurately. The computational domain extends from x = -0.015 to 9.0 
inches in the axial direction in this computation. Calculations were performed using a grid (32 x [127 x 256 \ 
1 million) depending on the size of the domain. Due to the very fine grid requirement near the nose, the allowable 
time step is very small and the computations become very expensive to simulate the unsteady computations in the 
entire domain. 

The acoustic field that impinges on the outer boundary is taken to be in the following form. 

p' = Real {pe‘ v±ie - y - imt } 

Here a ac , £ ac are the x, y wave numbers, respectively, of the acoustic wave and oo is the corresponding frequency 
of the acoustic disturbance. The incident angle 0 of the acoustic wave is defined as 

6 = tan -1 — 

a ac 

and in this paper computations are performed for zero incidence angle. 
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V. Results 


Computations are performed for hypersonic flows over cone and wedge with blunt leading edge, R n =0.001 in. 
Table 1 gives the flow parameters in the free stream and Fig. 1 shows the schematic diagram of the computational 
set up. First the mean flow field over a 5-degree straight cone with blunt nose tip is computed. After the mean flow 
is computed, acoustic disturbances are introduced at the outer boundary. The conditions at the boundary layer edge 
are given in Table 2. 


Table 1. Flow parameters for Horvath’s 15 wind tunnel model. 


Free stream 



Mach number 

M 00 = 

6.0 

Reynolds number 

Reoo = 

7.8xl0 6 / ft 

Density 

poo = 

7.059x1 (T 3 Ibm / ft' 

Velocity 

u 00 = 

3140.21 ft/s 

Reservoir Pressure 

Po = 

475 psi 

Reservoir Temperature 

T„ = 

475 °F 

Wall temperature 

Adiabatic condition 

Prandtl number 

Pr = 

0.70 

Ratio of specific heats 

y = 

1.4 

Length scale (x 0 =0.5 inch) 

“Vi 

8 

o* 

II 

7.30882xl0" 5 ft 


Table 2. Conditions at the edge of the boundary layer. 


(Var.) edee / (Var.)co 

Cone 

Wedge 

Mach Number 

0.932 

0.885 

Pressure Ratio 

1.560 

2.069 

Density Ratio 

1.372 

1.663 

Temperature Ratio 

1.137 

1.244 


Validation of Axisymmetric Solver. 

To validate the axisymmetric solver a series of comparison with the data available in the literature had been 
done. Flow field around a hemisphere-cone has been computed and results have compared with Mair’s 19 experiment. 
Table 3 gives the flow parameters and Fig. 2 shows the comparison of the flow field. 


Table 3. Flow parameters for Mair’s 19 wind tunnel model. 


Free stream 



Mach number 

M 00 = 

1.96 

Reynolds number 

Reoo = 

1.3xl0 7 / m 

Density 

Poo = 

0.2922 kg / m 3 

Velocity 

Uoo = 

511.76 m/ s 

Reservoir Pressure 

Po = 

14230.9 N / m 2 

Reservoir Temperature 

T 0 = 

169.65 °K 

Wall temperature 

Adiabatic condition 

Prandtl number 

Pr = 

0.70 

Ratio of specific heats 

y = 

1.4 

Length scale (x 0 =0.5m) 

ii 

b 8 

8 

6.5xl0" 3 m 


Ambrosio and Wortman 20 used experimental results to develop correlations for shock stand off distances as a 
function of Mach number. Figure 3 shows the comparison of shock stand off distances between the simulation and 
the predicted values from the correlation and Fig. 4 shows the comparison of shock shape at Mach 6.0 with 
Billing’s 21 correlation for cylinder-wedge and sphere-cone. It is seen the numerical simulation agrees with the 
correlations very well. 

Horvath et. al. 15 investigated the effect of bluntness on the transition onset at M=6.0 at different Reynolds 
numbers for a 5 degrees cone. They verified that the bluntness delays boundary layer transition in the hypersonic 
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flows. Computations were performed at three different nose radiuses as in the experiment and the wall to total 
temperature ratios are compared with the experimental results in Fig. 5. The figure shows a perfect match between 
the experimental and the simulation results. In the experiment, nose radius for the sharp cone is 0.0001 -inches and to 
perform a simulation with this small bluntness need very fine grid in the leading edge region and very small CFL 
number. To save the computational effort and time, a nose radius of Rn=0.001 -inches is chosen in this simulation. 
This explains the small shift in the temperature ratio between the computation and the experiment. 

Linear Stability 

Linear stability results for the similarity boundary layer over a cone and a wedge for the inviscid conditions at 
the surfaces are presented in Fig. 6. Figures 6(a) and (b) show the neutral stability diagram in (Re, F) plane for two- 
dimensional disturbances. The figure shows the first and second mode unstable regions and the variation of the wave 
number with the Reynolds number. Figures 6(c) and (d) show the N-Factor curves for different frequencies. Here the 
variables are non-dimensionalized by the variables at the edge of boundary layer. To obtain the variables non- 
dimens ionalized by the free stream values as given in Table 1, the variables in this section should be multiplied by 
the appropriate factors from Table 2. The frequency variable F has to be multiplied by 1.174 to obtain the values in 
terms of free stream values. 

The neutral stability curve clearly shows the unstable first and the second mode regions for the boundary layers 
over the cone and wedge at a free stream Mach number 6.0. The first mode and the second mode neutral stability 
curves merge at a Reynolds number of Re=1600 for the cone and they do not merge for the wedge case for these 
parameters. The most amplified frequencies are higher for the cone boundary layers than that for the wedge. The 
transition Reynolds numbers based on an N-Factor of 9 are about 3500 for the cone and about 5200 for the wedge. 
The most amplified frequencies are 0.85x1 0" 4 and 0.34x1 0" 4 for the cone and the wedge respectively. 


Acoustic Waves 

The linearized Euler equations in a uniform mean flow are: 




dt 


dx 
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The solution of this system can be written as 
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Here the pressure p is in the form 
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A Cl ac x+£ acy + Pac z - (0t ) 


p = p e 

i r amp 

The dispersion relation among the wave numbers CL ac , £ ac , j5 ac and the frequency CD is given by 

0 - <0 / = (<*J + + Pj W 


For acoustic disturbances with zero sweep (i.e. fi ac = 0 ), and zero incident angle (i.e. 6 = tan — 1 — = 0 ) the 

<*ac 

x-wave number a can be expressed as 


C O 

The plus sign corresponds to the fast moving wave and the minus sign corresponds to the slow moving wave. 

TT 03 

The corresponding phase speeds are C = U 0 ± a Q . The wave number of the fast moving wave is CL ac < 

U 0 +a 0 

CO 4 

and for the slow moving wave <2 > . For M = 6.0 and F = 1.2x10 , wave numbers of the fast and slow 

U o~ a o 

moving waves are a ac < 0.0586 anda ac > 0.0821 . Table 4 gives the wave number and the wave length for the slow, 
the fast frequencies. 


Table 4. Wave number and wave length for F=1.2xl0" 4 . 


Freq x 10" 4 



A - in 

1.20 

[467.79 kHz] 

Slow 

0.0821 

0.0671 

Fast 

0.0586 

0.0940 


Mean Flow for Cone and Wedge 

Figure 7 illustrates the mean flow data for the cone and the wedge computed using the WENO code. The figures 
on the left show the contours for the cone and the figures on the right show the results for the wedge. Figures 7 (a) 
and (b) show the mean flow density contours for the entire domain extracted from the Navier- Stokes computations. 
As expected the bow shock for the cone is narrower than that for the wedge because of the relieving effect of 
axisymmetry. Figures 7 (c), (d) show the density contours and (e), (f) show the Mach contours near the leading edge 
for the cone and the wedge respectively. The bow shock for the cone is located at 6 CO ne = 2.4xl0" 4 inch upstream of 
the leading edge and for the wedge it is at 6 we dge = 5.6xl0" 4 inch. Beyond the expansion fan the shock angles 
approach inviscid shock angle of 10.6 degrees for the cone and 13.1 degrees for the wedge. Figures 7(g) and (h) 
show the streamline patterns. The density profiles at x = 0.01, 0.03, 0.05, 8.94 inches. ( ^/Re x = 87, 151, 196, 2615) 

are plotted in Figs 8(a) and (b) in similarity coordinate for the cone and the wedge. The compressible Blasius 
similarity profile is also included for comparison. The boundary layer profiles approach to the similarity solutions 
close to x=0.1-inch for the cone. However, the density profile for the wedge at the end of the computational domain 
x = 9-in still did not approach the similarity results. This implies that the bluntness effects are stronger and persist 
for a longer distance for flow over wedges compared to flow over cones. Figures 8(c) and (d) depicts the same mean 
density profiles in y coordinate. Figure 9 (a) shows the mean pressure distribution along the surface for the cone and 

the wedge and Fig 9(b) shows the variation of the edge Mach number. At x = 8.9 -in M e = 5.57 and T =129.61 °R 

for cone and M e = 5.29 and T =141.79 °R for wedge. It can be seen that very close to the leading edge, there exists 

a strong shock and that the associated compression is followed by an expansion over the leading edge and the shock 
approaches the inviscid solution for the cone and the wedge. 
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Interaction of Slow and Fast Acoustic Waves with the Cone and Wedge Boundary Layer 

After the mean flow is computed two dimensional slow and fast acoustic disturbances are separately introduced 
at the outer boundary computational boundary and the time accurate simulations are performed. Unsteady 
simulations are performed for the frequencies F=1.2 xlO" 4 and 1.4 x 10" 4 for the cone case and for the frequency 
F=0.85 xlO" 4 for the wedge. These frequencies give maximum amplifications within the computational domain of 
x=9 inches. For the free stream parameters considered in these computations Table 1, these frequencies correspond 
to 331, 467 and 546 kHz respectively. To remain in the linear regime, the amplitude of the forcing freestream 

acoustic waves is given a small value ofp ac /p oo =2xl0 5 . Even with this small amplitude, nonlinearity starts to 
develop near the computational domain for the frequency F=1.2 xlO" 4 . 

Figure 10 shows the results for the evolution of unsteady fluctuations obtained from the simulations for the slow 
wave at a fixed time for the cone case with F=1.2 xlO" 4 . Figure 10(a) shows the contours of the density fluctuations 
in the entire domain and Fig. 10(b) depicts the results inside the boundary layer. In Figure 10(b) the surface of the 
cone is plotted along the x-axis to show the expanded view of the boundary layer clearly. The perturbation field can 
be divided into four regions. One region is the area outside the shock where the acoustic waves propagate uniformly. 
The second region is the shock layer across which the acoustic waves are transmitted. The third region is the area 
between the shock and the boundary layer. This region consists of transmitted external acoustic field and the 
disturbances that are radiated from the boundary layer. The fourth region is the boundary layer where the boundary 
layer disturbances evolve. The figures show that the disturbances inside the boundary layer are generated near the 
nose region. Figure 10(b) also shows the evolution of the first mode up to x~4.5 inches and the gradual 
transformation of the first mode to the second mode in downstream. Another interesting observation is that the 
region between the boundary layer and the shock layer is quieter compared to the acoustic waves outside the shock 
layer. This implies that the acoustic waves are weakly transmitted through the shock. This was also observed in the 
flat plate simulation (of Malik and Balakumar 15 ) where as the acoustic wave incidence angle is increased 
disturbances become quieter in the windward side. 

Figure 11 shows the evolution of the wall pressure fluctuations for the cone and the wedge case. Figures 11(a) 
and (b) show the pressure fluctuations induced by the slow and the fast acoustic waves for the frequency F=1.2 xlO" 4 
for the cone case and Fig 1 1(c) shows the results induced by the slow acoustic wave for the frequency F=1.4 xlO" 4 . 
Figure 1 1(d) shows the pressure fluctuation induced by the slow wave for the frequency F=0.85xl0" 4 for the wedge 
case. It is to be noted that different scales are used in Fig 11 due to the different amplification in the slow and the 
fast mode cases. For the frequency F=1.2 xlO" 4 , the maximum amplitude in the slow mode case is about 0.16 and is 
about 0.0025 in the fast mode case. This implies that the slow mode is more efficient; by about 67 times, in 
generating the instability wave inside the boundary layer compared to the fast wave. This agrees qualitatively with 
other simulations and analysis. One other observation is the amplification of the first mode near the leading edge 
region. The parallel linear computations revealed that the first mode is stable up to x~3 inch for this frequency. 
However, the simulation shows that the first mode disturbances are growing starting from the leading edge. Hence, 
the non-parallel effects are stronger in the cone case compared to the flat plate case and this yields a higher 
amplification ratio (about 67) between the induced flow field by the slow and the fast modes in the cone case 
compared to the flat plate case (about 20). Figure 1 1(d) shows that the amplification of the disturbances are small for 
the wedge case for this frequency F=0.85xl0" 4 . The maximum amplitude attained is about 2x1 0" 4 compared to 0.16 
in the cone case. The first mode region in the wedge case is stable contrast to the cone case. 

Figure 12 shows the amplitude of the pressure fluctuations along the wall in a log scale. This figure also includes 
the results from the parabolized stability equations (PSE) computations obtained for the same mean boundary layer 
profiles. The growth of the disturbances agrees very well with the PSE results. The figures clearly show the initial 
generation and the eventual exponential growth of the instability waves inside the boundary layer. The slow wave 
whose wavelength is closer to the wavelength of the instability wave transform into instability wave smoothly. The 
fast mode as was in the flat plate case initially generated the instability mode corresponding to the fast acoustic wave 
and then switches to the unstable second mode close to x = 4.5 inches. As was discussed earlier, the figure 12(a) 
shows that the first mode is growing due to the non-parallel effect. Due to the growth of the first mode starting from 
the nose region, it is difficult to define a receptivity coefficient in the cone case. We selected the amplitude near the 
nose region x~1.0 to compute the receptivity coefficient defined by the ratio between the initial amplitude of the 
pressure fluctuations at the wall near x~1.0 and the free stream acoustic pressure can be evaluated. For the wedge 
case, the first mode is decaying and a well defined neutral point is discerned. Hence, the receptivity coefficient is 
defined based on this amplitude. Table 5 gives the respective receptivity coefficients for the different cases. The 
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receptivity coefficients for the cone cases are about 4.569 and 0.068 for the slow and the fast acoustic modes and it 
is about 0.77 for the slow acoustic mode for the wedge case. 


c 


wall ) n 


recpt,p wall 


Table 5. Variation of receptivity coefficient for different X locations. 



Cone 

Cone 

Wedge 

F(xl0“ 4 ) 

1.2 

1.4 

0.85 

x(in) 

1.0 

1.5 

1.0 

1.5 

Neutral Point. 
(x=3.6) 

r 

recpt ,p wall ,Slow 

3.826 

4.569 

3.91 

4.89 

0.77 

r 

recpt, p wall .Fast 

0.057 

0.068 




r 

recpt, p wall ,Slow 

r 

recpt, p wall ,Fast 

67.12 

67.19 





Figure 13 shows the contours of the density fluctuations inside the boundary layer at three different streamwise 
locations to illustrate the structure and the evolution of the instability waves inside the boundary layer. The contours 
show that the disturbances are concentrated near the edge of the boundary layer and in downstream the disturbances 
exhibit the classical ‘rope’ like structures associated with the second mode. 

Discussion and Conclusion 

The receptivity and the stability of hypersonic boundary layers due to the interaction of two-dimensional slow 
and fast acoustic waves with a blunt 5 -degree cone and a blunt wedge are numerically investigated at a free stream 
Mach number of 6.0 and at a Reynolds number of 7.8x1 0 6 /ft. Both steady and unsteady solutions are obtained by 
solving compressible Navier-Stokes equations in cylindrical coordinates using the 5 th order accurate weighted 
essentially non-oscillatory (WENO) scheme for space discretization and using a third-order total- variation- 
diminishing (TVD) Runge-Kutta scheme for time integration. 

The shock stand off distances is computed for different bluntness at different Reynolds numbers and Mach 
numbers and is compared with experimental correlation results. The comparison showed excellent agreement 
between the simulation and the experimental results. The mean flow data for the cone and the wedge showed that 
the bluntness effect is stronger and persists longer for the wedge flows compared to the cone flows. 

The unsteady simulations showed that the instability waves are generated very close to the leading edge region. 
The simulations for the cone showed that the first mode starts to grow starting from the leading edge due to the 
nonparallel effects before they grow strongly due to the unstable second mode. In the wedge case, the first mode 
disturbances decay first before they start to grow due to the second mode. The receptivity coefficient of the 
instability waves generated by the slow acoustic wave is about 4 times the amplitude of the free stream acoustic 
wave. It is also found that the amplitude of the instability waves generated by the slow acoustic waves is about 67 
times larger than that for the case of fast acoustic waves. Therefore, forcing by slow acoustic wave is much more 
relevant in the transition process involved in hypersonic boundary layers. The receptivity coefficient in the wedge 
case is about 0.8 which is about 5 times smaller than that in the cone case. This is due to the initial growth of the 
first mode in the cone case compared to the decay in the wedge case and also due to the strong stabilization effect of 
the bluntness in the wedge case. The receptivity coefficient for a flat plate boundary layer with a bluntness of 0.0001 
inches 14 at a free stream Mach number of 4.5 is about 9 and the slow mode is about 20 times more efficient than the 
fast mode in generating the instability waves. This shows that the slow mode is much more efficient in flows over 
axisymmetric bodies than in two-dimensional flows. 
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Figure 1. Schematic diagram of the computational model. Parallel computation done with 32 block 
(32x[127x256]~lmillion points) 




Figure 2. (a) Comparison of Mach contours with Mair’s experiment 19 (b) Velocity vectors near stagnation 
point, (c) Streamline pattern colored by Mach contours near leading edge. 
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Figure 3. Comparison of computed and correlated 20 bow shock standing distance for sphere-cone and 
cylinder-wedge. 
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Figure 4. Comparison of computed and correlated 21 bow shock shape for sphere-cone and cylinder-wedge. White line represents 


the correlated bow shock shape. 
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Figure 5. Comparison of computed wall to total temperature ratio with conventional wind tunnel data 15 for 
different nose radii of 5 degree straight cone, (a) Rn=0.0625-in, (b) Rn=0.03 125-in, (c) Rn=0.001-in 

(d) Schematic diagram of the experimental model 15 . M = 6.0, Re = 7.8xl0 6 / ft,T = 475 °F,P q = 475 psi 
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Figure 6. Neutral stability diagram and variation of wave number with Reynolds number, (a) Cone, (b) 
Wedge (c) N-Factor curves for increasing frequency for cone and (d) for wedge. 
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results are on the right, (a), (b) Density contours in whole domain (c), (d) density contours, (e), (f) 
Mach contours, (g), (h) streamlines patterns near leading edge. 
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wedge respectively, (c), (d) density profiles versus Y coordinate for cone and wedge. 
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Figure 9. (a) Comparison of mean flow pressure distribution along the surface for wedge and cone, (b) 

Comparison of Mach number distribution at the edge of the boundary layer for wedge and cone. 
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X (in) x (in) 

Figure 10. Contours of the unsteady density fluctuations due to the interaction of slow acoustic wave with a blunt 
cone: F=1.2xl0" 4 . (a) Whole domain, (b) Fluctuations in the boundary layer. Computational domain 
rotated 5 degrees clockwise to show density fluctuations clearly. 



Figure 11. Wall pressure fluctuations generated by slow and fast acoustic modes, (a) F=1.2xl0" 4 Slow Mode- 
Cone, (b) F=1.2xl0' 4 Fast Mode-Cone, (c) F=1.4xl0‘ 4 Slow Mode - Cone, (d) F=0.85xl0‘ 4 Slow Mode-Wedge. 
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Figure 12. Amplitude of the pressure fluctuations at the wall and comparison with the PSE. (a) F=1.2xl0‘ 4 
Slow Mode, (b) F=1.2xl0‘ 4 Fast Mode, (c) F=1.4xl0' 4 Slow Mode, (d) F=0.85xl0‘ 4 Slow Mode. 
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Figure 13. Contours of unsteady density fluctuations inside the boundary layer. F=1.2xlO‘ 4 Slow Mode-Cone. R n =0.001 in. 
Computational domain rotated 5 degrees clockwise to show fluctuations in a better way. 
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